Isoscalar compression modes in relativistic random phase 

approximation 



Zhong-yu Ma^* , Nguyen Van Giai^, A. Wandelt^, D. Vretenar^^ and P. Ring^ 

1. China Institute of Atomic Energy, Beijing, P.R. of China 

2. Institut de Physique Nucleaire, IN2P3-CNRS, F-91406 Orsay Cedex, France 

3. Physics Department, Technical University Munich, D-85748 Garching, Germany 

(22 March 2000) 

O" 

O 

O. 

^ ■ Abstract 

A fully consistent relativistic RPA calculation is performed for the monopole 
^^ , and dipole compression modes in nuclei. The emphasis is put on the effects of 

Dirac sea states which are generally neglected in relativistic RPA calculations. 
■^ ' It is found that these effects can be quite important for the isoscalar monopole 

—1. ■ mode. The main contributions from the pairs of Fermi to Dirac sea states 

IQ ■ are through the exchange of the scalar meson, while the vector mesons play 

a negligible role. Numerical results of relativistic RPA are checked with the 
constrained relativistic mean field model in the monopole case. A good agree- 
ment beteween monopole energies calculated in RRPA and in time-dependent 
relativistic mean field approach is achieved. For the monopole compression 
mode, a comparison of experimental and calculated energies gives a value of 
250 ~ 270 MeV for the nuclear matter incompressibility. A large discrepancy 
^ I remains between theory and experiment in the case of the dipole compression 

H ' mode. 
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I. INTRODUCTION 

The linear response of a system to a weak external field can be calculated in the ran- 
dom phase approximation (RPA) which is one of the main theoretical approaches to study 
the nuclear giant resonances in a microscopic way |jl],H. In the non-relativistic approach 
the equivalence between the RPA and the small amplitude limit of the time-dependent 
Hartree-Fock approximation is well-known. Numerous RPA calculations based on Hartree- 
Fock single-particle spectra have been performed using effective nucleon-nucleon interactions 
such as Skyrme and Gogny forces. These calculations are very successful in describing the 
collective excitations and giant resonances in various nuclei in a framework which is consis- 
tent in the sense that both mean field and residual interaction come from the same effective 
force. 

In recent years the relativistic mean-field (RMF) theory with non-linear effective La- 
grangians has received considerable attention. The RMF has met great success in describ- 
ing the bulk properties of nuclei. One of the advantages of the RMF approach is that the 
same effective Lagrangian with only a few free parameters could predict the ground state 
properties of nuclei covering the whole periodic table and also exotic nuclei from the pro- 
ton drip line to the neutron drip line. The success of the RMF encouraged to investigate 
the extensions of the RMF approach to study dynamical processes, especially the collective 
giant resonances. The linear response of a system to an external field can be calculated in 
the relativistic RPA (RRPA) which is again a small amplitude limit of the time-dependent 
RMF (TDRMF) theory ||. 

In the relativistic approach, it is known from nuclear matter studies that self-consistency 
is extremely important. Self-consistent wave functions where the lower components in the 
Dirac spinor are enhanced by an effective mass must be adopted in the Dirac-Brueckner 
Hartree-Fock calculations |Q. Two questions of consistency appear in the RRPA. First, the 
particle-hole interaction must be determined from the same Lagrangian used in the RMF. 
The meson nonlinear self-interaction terms in the effective Lagrangian bring in density- 
dependent interactions which are essential for a quantitative description of the nuclear 
properties, especially the nuclear matter compression modulus. The meson propagators 
with nonlinear terms have been worked out in Refs. [|5|,0 . The second question of consistency 
is that the RRPA must be consistent with the RMF. It is assumed in the RMF that only 
positive energy states contribute to the nucleon self-energies (no-sea approximation). A com- 
plete relativistic Hartree approximation would have to include the contributions from Dirac 
states, which are divergent and a renormalization procedure is required ^. Unfortunately, 
a proper treatment of the relativistic Hartree approach in finite nuclei is too complicated 
and one has to work in the no-sea approximation. Nevertheless, a consistent RRPA built on 
an RMF ground state must deal in principle with the Dirac states. This is because RRPA 
is a small amplitude limit of TDRMF and the density change 6p{t) at time t has necessarily 
matrix elements not only between positive energy particle-hole pairs but also between pairs 
formed from Dirac sea states and occupied Fermi sea states of the static solution. Thus, 
an RRPA calculation consistent with the no-sea approximation should include both positive 
energy and negative energy single-particle states. Although it has been known before p| 
that the inclusion of the Dirac states ensures current conservation and decoupling of the 
spurious state, it was thought that the contribution to the RRPA energies and strengths 



from the Dirac states would be negligible. Therefore, the observables of interest (excitation 
energies and transition densities) would remain the same with and without including the 
Dirac states. Most previous works in studying the response functions of nuclei included the 
positive energy particle-hole pairs only. With this prescription, discrepancies between the 
RRPA and the TDRMF appeared, especially in the isoscalar monopole modes, which could 
not be understood for some time P,p!0|. The main purpose of this paper is to clarify this 
puzzle and to examine the effects of the Dirac states on the calculation of compression modes 
in the RRPA approach. We analyze the origin of the contributions from Dirac states in the 
study of the isoscalar nuclear giant resonances. To verify that the calculation of RRPA with 
Dirac states is the correct and consistent procedure we compare the inverse energy-weighted 
sum rules calculated in RRPA and in constrained RMF calculations. A good agreement 
can be obtained only when Dirac states are included. Results concerning isoscalar giant 
monopole resonances (ISGMR) and isoscalar giant dipole resonances (ISGDR) in nuclei are 
presented. 

The outline of this paper is as follows. The formalism of the relativistic linear response is 
outlined in Sec. II where we investigate the implications of including Dirac states in the RRPA 
equations. In Sec. Ill the constrained RMF method is used to calculate isoscalar monopole 
polarizabilities which are compared to those deduced from RRPA. Results of calculations 
of ISGMR and ISGDR in medium and heavy nuclei are presented in Sec. IV. Finally, Sec.V 
contains a summary of the main results. 

II. DIRAC STATES IN THE RRPA 

We start from the Hartree single-particle Green's function which is defined by: 

GH{x,y) = -t{0\T{^H{x)^Hiy)m , (1) 

where x = (t, r), "^nix) is the Hartree field operator and |0) is the uncorrelated Hartree 
ground state. Then, the Dyson equation for the nucleon propagator is 

Gh{x, y) = Go{x, y) + d'^xi(fx2Go{x, xi)J:^{xi,X2)Gh{x2, y) , (2) 

where the nucleon self-energy due to meson-nucleon couplings can be calculated in the 
Hartree approximation including only the tadpole Feynman diagram: 

S^(x, y) = S'^{x - y){-igl / d^XiD„{x - xi)Tr[GH{xi, x^)] 

-t^,g', J d'x.D^'^ix - x,)Tr[j,GH{xuxt)]} 
= S\x - y)[Es{x) - ^,E^,{x)] , (3) 

where x~^ refers to the time-ordering prescription and D^-, D^^ are the meson propagators of 
a and u, respectively. For the convenience of notations, we write explicitly only the a and 
uj mesons in this section but the actual calculations of the next sections contain also the p 
meson and photon fields. 

The Hartree Green's function in a standard spectral representation is expressed as: 



GH(ri,r2,E) = 2^ +}^ (4) 

+ E Mri)Mr2) (-^ - p \. ) 

h% \E-Sh-iV E-Sh + ir]) 

— Gh + Gh ! 

where we denote by a all positive energy states (e^ > 0), a all negative energy states 
{Sa < 0) and h the states in the Fermi sea. The unoccupied positive energy states will be 
denoted by p in the following expressions. The superscripts F and D refer to the Feynman 
and density-dependent parts of the Green's function, respectively. When one calculates the 
nucleon self-energy of Eq.(3) with this Green's function Eq.(4), all occupied states a and h 
will contribute. The contributions from occupied Dirac states are divergent and therefore, a 
renormalization procedure is necessary. This corresponds to the so-called relativistic Hartree 
approximation. In contrast, the RMF approach which is the framework we adopt here 
assumes no contributions to the nucleon self-energies from the Dirac states. This no-sea 
approximation corresponds to a modified single-particle Green's function, where the poles 
Sa + irj for negative energy states are shifted to Sa — it] in Eq.(4) [^. This gives a modified 
Green's function G^ built on a complete set of single-particle states and it satisfies the 
Dyson equation Eq.(^). A consistent RRPA built on the RMF ground state has to start 
from this modified single-particle Green's function. 

The RRPA method used in this paper is described in detail in Refs. |lTl|j5[| with just the 
replacement of Gh by Gh- The linear response of a system to an external field is given by 
the imaginary part of the retarded polarization operator: 

R{Q, Q- k, E) = -/mn^(g, Q; k, k; E) , (5) 

where Q is a one-body operator represented by a 4x4 matrix. The polarization operator 
can be obtained by solving the Bethe-Salpeter equation: 

n(P, Q; k, k', E) = no(P, g; k, k', E) (6) 

-T.9i J d'h£k2llo{P,r;kM,E)D,{k,-k2,E)U{T,,Q;k2,k',E) . 

i 

In this equation the index i runs over a, u and p mesons, Qi and Di are the corresponding 
coupling constants and meson propagators. The meson propagators for non-linear models 
are non-local in momentum space and they have to be calculated numerically. The detailed 
expressions of the -Dj(ki — k2, E) can be found in Refs. [^^. The Fj's are the 4x4 matrices 
1, 7^ and 7^r for i = a,uj and p, respectively. Notice that the space-like parts of vector 
mesons play no role in the static ground states, but they provide attractive contributions 
to the particle-hole residual interaction. The kernel of the Bethe-Salpeter equation has a 
simple form in the RMF approximation: 

no(P,Q;xi,X2) =2 {0\T[^Hix,)P^Hixi)^Hix2)Q^Hix2)]\0) (7) 

= i Tr[PGH{xi,X2)QGH{x2,xi)] . 

In the spectral representation, the unperturbed polarization operator has the following 
retarded form: 



HoV, Q; k, k', E) = Y. 

h, a=p,a 



[^h\P\^a){^a\Q\^h) ^ {^a\P\^h){^h\Q\^a) 



E -{sa- Eh) +iv E + {ea- Eh) + iv 



(8) 



which shows that the unperturbed polarization includes not only particle-hole pairs but also 
pairs formed from the Dirac and Fermi states ah, when one substitutes Gh in Eq-(7). The 
terms with ah pairs would correspond to the Pauli- violating terms which appear if one works 
with Gh in Eq.(7) [|ri|]. There is no Pauli violation in the present approach because we work 
with the modified Green's function Gh- 

In the time-dependent RMF (TDRMF) approach the no-sea approximation is imposed 
at each time. In the small amplitude limit the density p{t) can be expanded as the density 
at t = and a small density change Sp{t) which can be calculated in a complete set of 
states given by the static solution of RMF. This complete set includes the positive as well 
as negative energy states. Thus, Sp has non-zero matrix elements Spph and Spah- Therefore, 
the particle-hole space of RRPA must be built with both positive and negative energy states. 
The equivalence between the RRPA including the Dirac states and the small amplitude limit 
of the TDRMF with no-sea approximation can be formally shown [Q . 

In Eq.(8), one finds that the particle- hole energies between the occupied positive energy 
states and the Dirac states are usually two orders of magnitude larger than those involving 
only positive energy states. Furthermore, the matrix elements of P and Q between positive 
and negative energy states should be small because the positive (negative) energy states 
have a small lower (upper) component in their Dirac spinors. It is thus found that the ah 
pairs make no visible contribution to the imaginary or to the real part of the unperturbed 
polarization in the energy region of the nuclear collective excitations. This is the reason why 
the Dirac states were usually not included in most previous works of the RRPA. However, 
we will see in the next sections that the correlated polarization operator can be affected by 
the Dirac states when the particle-hole interaction depends on the a meson, i.e., in the case 
of isoscalar excitations. 

III. CONSTRAINED RMF 

In the non-relativistic approach there is a complete equivalence between the static polar- 
izability calculated in RPA and that deduced from a constrained Hartree-Fock calculation, 
and this property can be used as a check of consistency of the calculations. The same 
property will be used here to compare the static polarizabilities of RRPA and constrained 
RMF and to investigate the importance of the inclusion of Dirac states in the RRPA. The 



constrained RMF method was used in Ref. [12| to evaluate the ISGMR energies within the 



Walecka's linear model. In this section we discuss the constrained RMF results obtained 
with effective lagrangians having different types of self-interactions in the a and u fields. 
In general, a constrained RMF equation with an external field XQ is described as: 

{H{X) + XQ)\X) = E{X)\X), (9) 

where H(X) is the RMF hamiltonian, which depends on the parameter A due to the self- 
consistency of the single-particle states, |A) andE{X) are the ground state and corresponding 
energy of H{X) + XQ, respectively. At A=0 the first order variation of {X\Q\X) with respect 



to A, i.e., the static polarizability is equal to the second order variation of E{X). The same 
static polarizabihty can be calculated as twice the inverse energy-weighted sum rule m_i of 
RRPA, 

d{X\Q\X) _ ^^Ex _ 
dX '^"° dX^ '^"° 

_o^ (0|Ql^n)(v&n|g|0) ^ ^^Q^ 

The second line of Eq.(^) shows that the polarizability can be deduced from the real part 
of n calculated at E=0. In the case of the ISGMR, the operator Q is simply 7or^ and the 
RRPA m_i sum rule can also be obtained by performing a spherically constrained RMF 
calculation. Note that, if one replaces in Eq.(^ the self-consistent hamiltonian H{X) by the 
fixed hamiltonian Hq of the unconstrained RMF solution, then Eq.(|lO|) gives the m_i sum 
rule of the unperturbed response corresponding to the polarization operator Ho. 

The effective Lagrangians we use here include the scalar meson a, vector meson u, 
isovector vector meson p and the photon field A^^. Self-interactions in the a and uj fields are 
introduced in order to simulate a density dependence and hence to lower the compression 
modulus of nuclear matter 1 11 311: 



C = ^Hiz^d, -Mm- g.cr - g^^uJ, - g^r'^^Pl ~ 67^,^(1 - Tg))* (11) 

where 

W^"" = d^u" - d'^uj" , (12) 

Ua = -mla'' + -g^a^ + -g^a\ U^ = -m>^a;^ + -Csiuf^u^f . (13) 

We briefly outline the procedure of the numerical calculations. The nucleon mean fields are 
calculated by solving the Dirac-Hartree equation in coordinate space self-consistently in the 
RMF approximation. All single-particle states (unoccupied and occupied) are calculated 
by diagonalizing the self-consistent mean fields in a harmonic oscillator basis with a proper 
choice of the oscillator length b and the number of shells (A^=12). An averaging parameter 
A=2MeV is introduced as the imaginary part of the energy to smooth out the response 
functions. Calculations are performed with different non-linear lagrangians which span a 
wide range of values of the nuclear incompressibility: NLl [Q, NLSH [|T3[, NL3 |]TB| and 



TMl |T^. They generally give good descriptions of the ground state properties. The linear 
Walecka model HS |T^ is also considered for comparison because of its high incompressibility. 
In the constrained RMF calculations we use a small step AA=0.5xlO~^fm~^ in order to 
calculate the derivatives with sufficient accuracy. 



The values of the inverse energy-weighted sum rules calculated in constrained RMF 
and in RRPA are compared in Table 1 for the case of ISGMR in ^''^Pb. We find a good 
agreement (within 1%) between the first two columns, which is a numerical check of the 
accuracy of the constrained RMF calculations. The values of m_i in the third column are 
calculated by integrating the inverse energy- weighted ISGMR strength (from to 60 MeV). 
The last column corresponds to one half of the real part of the response function at zero 
energy. The agreement between the last two columns simply indicates that the dispersion 
relation between the real and imaginary parts of 11 or Ho is well obeyed numerically. The 
comparison between the constrained RMF and the RRPA provides a severe check for the 
RRPA method used in this paper. It is found that the agreement is quite good (within 
1%) for the unperturbed cases while the difference becomes a few percent for the RRPA 
cases. We find that the discrepancies between RRPA and constrained RMF are mainly due 
to the Coulomb force. If one switches off the Coulomb force the agreement becomes perfect 
(less than 1% in all cases). The expansion of the wave functions on the harmonic oscillator 
basis may not be accurate enough to account for the long range effects of the Coulomb force 
especially for the Dirac states with an attractive Coulomb interaction. A more sophisticated 
method would be to calculate the Green's functions directly in coordinate space ||T9| . 

The discrepancies between the RRPA and the constrained RMF are dramatic if one 
does not take the Dirac states into account as one can see in row III of Table 1. The 
large contribution to the RRPA strengths from the Dirac states can be understood in the 
following way. In the relativistic approach the nucleon potential is obtained by a strong 
cancellation between the isoscalar scalar and vector potentials, and so do the particle-hole 
residual interactions in the isoscalar channel. Due to the 7-matrix structure the matrix 
elements of the residual interaction between particle-hole pairs through the exchange of 
scalar and vector mesons are largely cancelled to two orders of magnitude. However, this 
is not true for pairs built with Fermi and Dirac states because the large component of the 
wave function is the upper component (lower component) for a Fermi (Dirac) state. The 
orthogonality between positive and negative energy states with the same quantum numbers 
is mainly due to the cancellation between the upper and lower components. Again, due 
to the 7-matrix structure the matrix elements between Fermi and Dirac states through the 
exchange of vector mesons {uj or p) are small whereas those corresponding to a scalar meson 
are large. These scalar meson matrix elements are repulsive and they can overcome the large 
energy gap between Dirac and Fermi states. 
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FIG. 1. ISGMR strength distributions in ^*^^Pb calculated with NL3 paranietrization. The 
solid and long-dashed curves are the RRPA strengths with and without Dirac states, respectively. 
The short-dashed (dash-dotted) curve corresponds to calculations where only vector mesons (scalar 
mesons) are included in the couplings between Fermi and Dirac states. 

The large effects of Dirac states are illustrated in Fig.l. We separate the contributions 
from scalar and vector mesons for the ah pairs and we plot the monopole strengths in ^°^Pb, 
taking NL3 as an example. The solid and long-dashed curves are the RRPA strengths with 
and without Dirac states, respectively. They are markedly separated. The ISGMR strength 
is pushed strongly down without the Dirac state contributions. The short-dashed curve is 
the RRPA strength with the contributions of ah pairs coming from the vector meson only. 
It is more or less the same as the case without Dirac states. The dash-dotted curve is 
calculated with only the scalar meson contributions to the ah pairs. It is clearly seen that 
the contributions of the scalar meson in the ah pairs dominate whereas those due to the 
vector meson are negligible. 

It can be concluded that the inclusion of only the positive particle-hole pairs in the 
isoscalar modes, where the isoscalar mesons play an important role, would provide a too 
strong attraction and therefore, a too low isoscalar giant resonance energy. This will reflect 
in a too large polarizability as shown in Table 1. On the other hand, the contributions 
from Dirac states can be neglected in the case of isovector modes where the isovector- vector 
meson dominates. 



IV. MONOPOLE COMPRESSION MODE 



The value of the compression modulus Koo of nuclear matter is one of the most impor- 
tant issues of nuclear physics. The generally accepted procedure for determining K^o is to 
calculate isoscalar monopole energies, as well as other properties of finite nuclei using various 
microscopic models in order to select those which can describe correctly the data and to use 
these models for constructing the nuclear matter equation of state [^]. Here, we examine 
the predictions of the various effective lagrangians in the framework of RRPA. In this work 



we adopt the mean field approach and therefore the effects of 2p-2h couplings |2T[ are not 
included. 
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FIG. 2. ISGMR strength distributions in Pb calculated with different parameter sets. 

We are interested in establishing a link between the measured energies of compression 
modes in nuclei and K^o- The most favorable situation is met in heavy nuclei where the 
giant monopole resonance is less broad than in medium and light nuclei. In the latter case, 
the Landau fragmentation is quite strong due to the weakness of the particle-hole interaction 
in the isoscalar monopole channel. In order to have a comparison with non-relativistic RPA 
results and with TDRMF results we concentrate on ^osp]^ ^^^^ ^^^Sm where the experimental 
giant monopole energies are 14.2±0.3 and 15.4±0.3 MeV, respectively [^. The calculations 
are performed with 5 effective Lagrangians whose compression modulus span a wide range 
of values: NLl {K^ = 211 MeV), NL3 (i^oo = 272 MeV), TMl {Koc = 281 MeV), NLSH 
{Koo = 355 MeV), HS {Koo = 545 MeV). The inclusion of negative energy states in the 



RRPA calculations necessarily involves the choice of a cut-off energy in the Dirac sea. The 
bottom of the potential on the negative side is around -150 MeV, typically. We have checked 
that the results are stable if we move the cut-off between -300 MeV and -900 MeV. 

The RRPA monopole strengths with and without Dirac states are shown in Fig. 2 for 
different effective Lagrangians. The strengths are shifted down strongly to the lower energies 
when only the positive energy particle-hole pairs are included. In some parametrization, 
such as NLl, the strong attractive residual interaction produced by the positive particle- 
hole pairs would bring the system beyond a critical point, where the system is not stable. 
In the complete RRPA calculations with the Dirac states included the strengths are shifted 



to higher energies and the unstabilities disappear. In a recent work [^] such an effect due 
to the Dirac states has been shown analytically in the framework of nuclear matter and the 
Walecka's linear model. 
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FIG. 3. RRPA peaks and centroid energies of the ISGMR in ^osp]^ ^^^^ ^^^Sm as functions of 
Koo- TDRMF results for ^''^Pb[24] are also shown. The experimental values for ^osp^-, ^.^^^ ^^^Sm 
are E = 14.2±0.3 and 15.4±0.3 MeV, respectively [22]. 

In Fig. 3 we show the RRPA monopole energies in ^°^Pb and ^^"^Sm as functions of Koo- 
In each case, the difference between peak and centroid energies is generally small due to the 
fact that the strength is mostly concentrated in the giant resonance region. This difference 
would be larger in lighter nuclei. The ISGMR energies calculated in TDRMF [^ are also 
plotted in figure 3 in comparison with our RRPA results. Note that the present results 



differ from those of Ref. |^ where the space-like contributions of vector mesons were not 



10 



included in the particle-hole residual interaction. A good agreement of the ISGMR energies 
from the two approaches is obtained. The dependence on K^o is nearly linear, similarly to 
what was obtained in non-relativistic RPA calculations performed with Skyrme and Gogny 



interactions [^. From Fig. 3 it can be concluded that a value of K^o close to 250 ~ 270 



MeV would bring the RRPA results in good agreement with experiment. 

V. DIPOLE COMPRESSION MODE 

The isoscalar giant dipole compression mode, or squeezing mode is mainly built on 3/ia; 



excitations. It was already observed a long time ago p^ and it has been measured recently 



with increased accuracy by (a, a') scattering |2^ in several medium and heavy nuclei. We 



have performed RRPA calculations of this mode in ^°^Pb. An external field [^ 7o(^^~?7^)^io 
with ?7 = I < r^ > is used as a probe. The second term in the operator is introduced to 
remove a large amount of spurious components coming from the admixture of the center-of- 
mass state. The value of i] is determined by a condition of translational invariance. 

It is found that there is a large amount of strength in the Ihuo region around 10 MeV. 
The main peak (in terms of energy-weighted strength) is at about 26 MeV for the case of 
NL3. This peak moves up with increasing incompressibility. This general pattern of results, 
namely a sizable fraction of strength at lower energies and a giant resonance peak well above 



20 MeV is also found in non-relativistic RPA calculations with Skyrme interactions P9| , |21 
This is at variance with the experimental findings [^ where no strength is observed in the 
10 - 15 MeV region and where the ISGDR peak is close to 20 MeV in 2°^Pb. It should be 
noted, however, that in the continuum-RPA calculation of ref. [^ the low-lying strength 
could be eliminated by a particular subtraction procedure, but the ISGDR peak still remains 
above 20 MeV. 

In Fig. 4 are shown the calculated peak energies as a function of the compression modulus. 
The centroid energies calculated in the giant resonance region (20-40 MeV) are also displayed. 
The peak and centroid energies are generally quite close except for the NLSH case where 
the giant resonance is much Landau- fragmented. The calculated energies have a weaker 



i^oo-dependence than in the monopole case. The measured value |]28| of centroid energy is 
20.3 MeV. One can see that a model like NL3 which is doing rather well for the ISGMR 
energy is now off by 5 MeV in the ISGDR case. This disagreement between experiment and 
models for the ISGDR case is the next puzzle to be solved. 
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4. Peak and centroid energies (in the 20 - 40 MeV region) of the ISGDR in ^''^Pb as a 
of i^oo- The experimental centroid energy is E = 20.3 MeV. 



VI. SUMMARY AND CONCLUSION 

In this paper we have studied an RRPA approach built on the RMF in the framework of 
the relativistic baryon-meson field theory. The effective Lagrangians with non-linear meson 
self-interactions, which describe well the ground state properties of finite nuclei as well as 
the nuclear matter matter properties are adopted. The calculation is fully consistent in the 
sense that the single-particle basis and the particle-hole couplings are obtained from the 
same effective Lagrangian. In a no-sea approximation, a consistent treatment of the RRPA 
must include both the particle-hole pairs built with positive energy states and the ah pairs 
formed with the occupied Fermi sea and the empty Dirac sea states. 

The importance of the inclusion of the Dirac sea states is examined in detail in the case 
of the monopole compression mode. We have used the constrained RMF method with a 
r^ constraint as a check of RRPA calculations. With the inclusion of Dirac sea states the 
RRPA static polarizability, or the inverse energy-weighted sum rule, is in agreement with 
the constrained RMF result whereas the RRPA polarizability becomes much larger (it can 
even be negative for some effective Lagrangians) if the Dirac states are omitted. In general, 
the Dirac sea states must play an important role in the isoscalar modes. The reason is that 
the repulsive contributions of the ah pairs through exchange of the isoscalar-scalar meson 
dominate. Their effect is to shift the isoscalar giant resonances to higher energies. 
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The remaining small discrepancies between RRPA and constrained RMF polarizabilities 
when the Coulomb field is turned on seem to indicate that the harmonic oscillator basis 
used here may not be fully appropriate. A proper treatment of the Coulomb interaction 
in Dirac sea states for heavy nuclei is required to get accurate results because of the long 
range attraction for the negative energy states. Employing the Green's function technique 



in coordinate space and carrying out a continuum RRPA [^ will probably improve the 
accuracy of results. A good agreement of ISGMR energies calculated in RRPA and TDRMF 
is obtained when we include the Dirac state contributions in the calculation of RRPA. 

The comparison of the isoscalar monopole data for heavy nuclei, namely ^ospb and ^"^^Sm, 
with calculations performed with different parameter sets indicates that the compression 
modulus of nuclear matter is in the range of 250 ~ 270 MeV. These values are somewhat 



larger than those predicted by non-relativistic RPA models p6 |. 

The main open problem is now to understand the ISGDR energies. The discrepancies 
between experiment and predictions for this mode are the largest ever found for isoscalar 
giant resonances described in non-relativistic RPA or RRPA. It is unlikely that the coupling 



to 2p-2h configurations [^ can account for them 
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TABLES 

TABLE I. The ISGMR inverse energy-weighted sum rule in ^''^Pb. The rows with I, II, III 
refer to the results of unperturbed (Hartree), RRPA with Dirac states and without Dirac states, 
respectively. All values are in lO^fm^MeV"^. 
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